Seasonal Composite Landsat TM/ETM+ Images Using the Medoid (a Multi-Dimensional Median)

Multi-temporal satellite imagery can be composited over a season (or other time period) to produce imagery which is representative of that period, using techniques which will reduce contamination by cloud and other problems. For the purposes of vegetation monitoring, a commonly used technique is the Maximum NDVI Composite, used in conjunction with variety of other constraints. The current paper proposes an alternative based on the medoid (in reflectance space) over the time period (the medoid is a multi-dimensional analogue of the median), which is robust against extreme values, and appears to be better at producing imagery which is representative of the time period. For each pixel, the medoid is always selected from the available dates, so the result is always a single observation for that pixel, thus preserving relationships between bands. The method is applied to Landsat TM/ETM+ imagery to create seasonal reflectance images (four per year), with the aim being a regular time series of reflectance values which captures the variability at seasonal time scales. Analysis of the seasonal reflectance values suggests that resulting temporal image composites are more representative of the time series than the maximum NDVI seasonal composite.